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For the Casimir piston filled with an inhomogeneous medium, the Casimir energy is regularized 
and expressed with cylinder kernel coefficients by using the first-order perturbation theory. When 
the refraction index of the medium is smoothly inhomogeneous (i.e., derivatives of all orders exist), 
logarithmically cutoff-dependent term in Casimir energy is found. We show that in the piston model 
this term vanishes in the force and thus the Casimir force is always cutoff-independent, but this term 
will remain in the force in the half-space model and must be removed by additional regularization. 
We investigate the inhomogeneity of an exponentially decaying profile, and give the first-order 
corrections to both free Casimir energy and Casimir force. The present method can be extended 
to other inhomogeneous profiles. Our results should be useful for future relevant calculations and 
experimental studies. 


I. INTRODUCTION 

Casimir effect!^ is known as one of the direct man¬ 
ifestation of the vacuum zero-point energy in quantum 
physics. A mode-summation method can be used to pre¬ 
dict easily an attractive force between two electrically 
neutral, perfectly conductive half-spaces. This ghostly 
force arises from the quantum fluctuations, and thus 
is ubiquitous in the physical world. However, once 
predicted, the Casimir force in a variety of scenarios 
has witnessed many divergence problems during calcu¬ 
lations, due to the geometry of boundaries or topology 
of spaced—. These unresolved divergences usually show 
a logarithmic cutoff dependence and seem to be irremov¬ 
able (while other quartic or cubic cutoff-dependent terms 
in Casimir energy are already well understood as volume 
energy or surface energy, and thus are removable) 6,7 . 

When the inhonrogeneity of the medium rather than 
complicated geometry of objects in Casimir apparatus is 
considered, it is known that analytical description of the 
Casimir force (stress or force density) has already been 
obtained^, without explicit divergence problem. The 
results in both references have subtracted the “bulk con¬ 
tribution” and thus relate only to the scattering Green’s 
function. This, in our minds, regularizes the quartic 
diverging term (volume energy) in Casimir energy, but 
might not be a thorough regularization, as also men¬ 
tioned in Ref<£d£. In fact, divergences indeed occur when 
one tries to evaluate the numeric values of both Casimir 
energy and Casimir force, according to the obtained an¬ 
alytical description within an inhomogeneous medium^. 
Efforts have been made, since then, to understand the 
remaining divergence^, but the problem is still far from 
solved. 

We note, in the model of half-spaces which is used in 
Ref.®, cutoff-dependent terms in Casimir energy must be 
assigned physical meanings and thus removed manually 
by introducing corresponding additional regularizations. 


While in the model of Casimir piston, cutoff-dependent 
terms may vanish automatically in the Casimir force, due 
to the cancellation of contributions from the left and 
right cavities. We also note, a naive mode-summation 
approach!- 2 in the first-order perturbation theory (which 
is supposed to show the cutoff dependence), turns out 
to yield cutoff-independent result of Casimir force for in¬ 
homogeneous media. All of the above have led us to ex¬ 
pect that more useful information could be obtained from 
the Casimir piston model. Therefore, we adopt the pis¬ 
ton model here, following the mode-summation approach 
in the first-order perturbation theory as Refi!2(but in a 
more general form), to investigate the Casimir physics 
within inhomogeneous media. The purposes for do¬ 
ing so are threefold. First, we want to analyze the 
inhomogeneity-induced Casimir divergence with the heat 
kernel expansion, and expect to obtain some insights for 
additional regularizations for half-space model. Secondly, 
we want to show that the cutoff independence of Casimir 
force is true for various inhomogeneous profiles, instead 
of a particular case. Thirdly, we want to show how large 
the influence of the inhomogeneity on the total Casimir 
force is. As weak force measurements have been devel¬ 
oped and reached a quite precise level (within 1%)!^, 
and Casimir force between bodies in a liquid has also 
been measured ! 4 ' ^ ,, inhomogeneity-induced corrections 
may be useful for future experimental studies as the ex¬ 
perimental configuration becomes more and more com¬ 
plicated. 

In the present paper, we first derive general expres¬ 
sions of Casimir energy for inhomogeneous media in 
Sec. El The summations in expressions of Casimir energy 
are organized and re-expanded over the cutoff parame¬ 
ter in Sec. nn Then we prove the cutoff independence 
of Casimir force for smooth inhonrogeneity in Sec. m 
where we also see logarithmic divergence in the half-space 
model. We investigate the exponentially decaying inho¬ 
mogeneity in Sec. [V] and give first-order corrections to 
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FIG. 1. Casimir piston model^~ with four (up, bottom, front 
and back) hidden plates. The cyan curve represents an arbi¬ 
trary inhomogeneous profile. All plates are of perfect conduc¬ 
tivity. 


Casimir energy and Casimir force. Discussions and con¬ 
clusions are given in Sec. |VT]and Sec. IVII1 respectively. 


II. CASIMIR ENERGY OF PLATES WITHIN 
INHOMOGENEOUS MEDIA 

In the mode-summation technique, total zero-point en¬ 
ergy of the Casimir piston device, see, Fig.Q], is expressed 
as^ 

£■0=2 + £ —> (1) 

m,p,q, X 

where m,p,q are indexes for three wave numbers (m for 
the direction perpendicular to the plates), A is the index 
for polarization, notation C —> 1Z represents the counter¬ 
part for the right cavity, and also we have set h = c = 1 
which will be recovered later according to dimensional 
analysis. For simplicity, we use fc|| to account for {p, g}, 
t for {to, fc||} and J for {f, A}. Below we will focus on 
the left cavity and omit the notation C —> TZ in all equa¬ 
tions, just keeping in mind that the right counterpart 
should be added in the final step. We adopt the cutoff 
regularization, Eq. |Tj then becomes 



We omit the limit notation in the following for simplicity 
as well. 


When there is an inhomogeneous perturbation in re¬ 
fraction index n(x) = no[l+5a-f(x)], where 8a is a small 
perturbation value, the difference of regularized Casimir 
energy is then 

^ = ^E^ e_Ch,0j > ( 3 ) 

j 

with u) = -6aVjU}°j and Vj = (x°j\f (x)\x°j) ■ Xj is the 
J t h unperturbed eigen-wavefunction with being its 
eigen-frequency. We can find the expressions for the elec¬ 
tric fields Xj in RefJ^ (Eqs. (14) and (15) after correction 
in their erratum; we do not use E for electric field here to 
avoid confusion with Casimir energy). We have also con¬ 
sidered the perturbation of frequency in the exponent, 
and thus we get a factor of ■ £• in our Eq. [3] However 
it does not otherwise change our argument. Therefore, 
we can obtain 

5E = -^d r tY / Vjk 0 je~t k ° /n °. (4) 


Here kj = \J^\\ + ( m7r /a) 2 is the wave number for homo¬ 
geneous media defined as fcj = noWj/c, and we have set 
c = 1. The perturbation theory is justified as long as Vj 
is bounded and daVj <C 1. We note that only Vj varies 
for different polarizations. To sum up polarizations first, 
we obtain 




(2-8 m0 )F 0 -(2-S m0 ) 


(tott/o) 2 
A:| + (rmr/a) 2 ’ 


( 5 ) 


with the Fourier coefficient of perturbation profile f(x) 

„ 1 f“ 2rmrx . . 

Fra = - f(x) COS - dx, (6) 

a J o a 

where A = 1,2 represents two different polarizations. To 
obtain the counterparts of our Eqs. [5] and [G] for the right 
cavity we can substitute a with D — a and also change 
the integration range in Eq. [G]from (0,a) to ( a,D ). 

Substituting Eq. [5] into Eq. [4j we can split the Casimir 
energy into two parts for the convenience of calculation, 
the first part is 

SE, = (2 - 6 m0 )F 0 k?e-Z k °/ n ° 

2n 0 ^ 

= ^^£o^S(0 £ ( 2 - Smo) e~ m ^ an \ 

m 

( 7 ) 


and the second part is 


8a 


SE 2 = —^£(2-*, 


2 n, 

SaArig 
47r 


(TOT Tjaf_ o e -tk°/no 

mOy 7.2 , /_ / ..\ 2 A 


fcn + (rmr/a) 


cf E ( 2 - 6 ™o)F„ 


0 —m'K^/ano 


( 8 ) 
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where S = £ _2 (1 — £3f), m = 0,1, 2, • • ■, and A is the 
surface area of plates. To get the second equalities for 
both Eqs. [T] and [8] we have integrated over k\\, provided 
that the lateral dimension of the plates is much larger 
than the distance between plates, i.e., y/A >> D, so we 
can replace summation with the corresponding integral. 
The counterparts of Eqs. |3 and [8] for the right cavity can 
be obtained by substituting a with D — a. 

On the other hand, following Eq.[2]and integrating over 
k ||, we can also obtain the regularized Casinrir energy for 
the homogeneous case 

/t ^.2 _ 

E h = —^0 € H(£) £ (2 - 5 m o)e~ m7V ^ an °. (9) 

m 


while we have F^ 0 = f£ : = 0. The explicit expression 
is then 


rpC _ 
r m,d 


(-1) 1 1 d\ 2 i /a\ 2i ~ 1 /a\ 

T ( d + i-2»)! (2mjr) (s) U 

i=l 


d+l-2i 


(14) 

where [d/2] denotes the bare integer part. The counter¬ 
part for the right cavity can be obtained similarly 


[d/2] 


F n = V" ( ~ 1)a lrf! (2rmr)~ 2i 
m ’ d ^ {d+l-2i)\ [} 

i=1 



1 ( a \d+l-2i 

V D ) 

[ D > 


(15) 


Now, we can multiply Eqs. [3 [8] and [9] with a factor of 
he, and add to them with the counterparts of the right 
cavity, and then take the £ —> 0 limit. The sum of these 
three equations gives the total Casimir energy for inho¬ 
mogeneous case. Its derivative with respect to position a 
yields the Casimir force on the central plate. These three 
expressions are quite general for perturbation profiles as 
long as Vj is bounded and 5aVj <C 1. 


We note that they all have the form where 

<7i’s are some coefficients independent of to. Back to 
Eq. El we now know the summation is a combination 
of polylogarithm functions 

c a 9 [c^/2] 

sE 2 = E 9i • Li 2i (e~^ an °) . (16) 

i=l 


III. SUMMATION AND RE-EXPANSION 

The summations in Eqs. [3 and [3 can be expressed in 
terms of polylogarithm function 

E(2- 6mo)e- mni/ano = ( E 2e~ m ^/ an ° ) + (2 - d 0 o) 

m =0 \m =1 / 

= 2 Li 0 + 1. 

( 10 ) 

The summation in Eq. [8] depends on the particular form 
of the perturbation profile. However the to = 0 term 
vanishes obviously. Then we can rewrite it as 

SE 2 = S -^d! F m e~ m ^ /an \ (11) 

m—Oorl 


In general, these polylogarithnr functions can be ex¬ 
panded over £ (with £ —> 0) as 



where s = 2,4,6, •••, F[ s = l//i is the harmonic 

number with Hq = 0. 

Before going further to get the explicit expression of 
Casimir energy for inhomogeneous case, we inspect an¬ 
other set of profiles 


f{x) = e~ vx+A , (18) 


where to can run from 0 or 1 as needed. For the Taylor 
basis of order d 


where Re[r]\ > 0. For the left cavity we have (we now 
consider m running from 0) 


M*) = (§)", (12) 

where d = 0,1,2, • • •, according to the definition of F m 
(we now consider m running from 1), we can integrate by 
parts to get the following relation for d > 2 for the left 
cavity 

c ad ( a ^d-i , a n 2 d{d-l) c 

m ’ d D(2mir) 2 ' D l 2m7C D 2 m ' d ~ 21 

(13) 


F 


m,rj 



— e 


—rja 


l 

to- — iija/2n 


1 

to + irja/ 2-7T 


(19) 

The right counterpart could be obtained by substitution 
(1 — e~ va ) —> (e~ va — e~ vD ),a —>■ (D — a). Eq. [Til then 
becomes 




•«? 


27r 

4>{e" 


47r* 


i lr l a \ ~ 


^,1,—) 


2n ' 


( 20 ) 
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Here Lerch zeta function <f> can be expanded over £ (with 
£ —» 0) as 


<t>(e “"°,l,/3) = 


+ 


(JlLff)k 


. fe =0 


fc! 


f—JLL) k 

£ca 

k= 1 


k\ 


^(1) - V>(/3) - l°g (-) - !«g£ 

ano 


( 21 ) 


where ip is the digamma function and £ is the Hurwitz 
zeta function. We note /3 is not any negative integer and 
thus £(1 — k,/3),ip(/3) are always finite. Furthermore, 
since only /3 is complex, we have 0(e _7r ^' /an °, 1,/?*) = 
</>*(e -7r £' /ano , l,/3) and that is why Eq. [20] is always real. 


IV. CUTOFF INDEPENDENCE 

Up to now, we see the total Casimir energy for inhomo¬ 
geneous case is generally expressed as a Laurent-type-like 
series 

E = E h + 6E, + 6E 2 = J2 Cit + 51 Cl ^ lo §£- ( 22 ) 

i=— 4 2=0 

These coefficients C,; and Cli are well-studied heat ker¬ 
nel (or more precisely cylinder kernel here) coefficients 
in cutoff regularizatio n 6 i 7 i 16 . They depend only on the 
geometry property and the boundary condition of the 
system under consideration. The divergent terms under 
limit £ —► 0 are usually assigned to the self-energy of vol¬ 
ume or surface and so on to renormalize the theory. Here, 
in our case, we should check these divergent terms, mak¬ 
ing sure that their coefficients are independent of position 
a so that these divergences would not go into the Casimir 
force, as we expect the observable—the Casimir force— 
to be finite. The constant term Co is the free Casimir 
energy, and its derivative with respect to position a is 
the Casimir force we want to calculate. Once these di¬ 
vergent terms are assured to be independent of position 
a, the final result of Casimir force should be identical to 
the result from Ref; 17 . 

We now prove the cutoff-independence of Casimir force 
for any smoothly inhomogeneous perturbation profile. 
We know, any smooth function in domain (0, D) can be 
expanded via the Taylor bases given in Eq. Q21 There¬ 
fore, equivalently what we need is to prove the cutoff- 
independence of Casimir force for basis of any order d. 
We proceed in the following way. 

First, the cutoff property of Eq. 0 is determined by 
Fo<9{2(£)(2L*o + 1) and the cutoff property of Eq. [5] is 
determined by <9|S(£)(2 Li 0 + 1) according to Eq. [TUI We 
note the operators acting on the polylogarithm function 
generally decrease the power of £ by three orders, and 
Li 0 = X)i=-i hC- Thus we have to check i = — 1 to 3 
to see the cutoff property. According to the expansion of 


polylogarithm function, Eq. 1171 we have 

,c ,-i _ an 0 , n ^-i _ (D - a)n 0 

l -is — —— — ■ 

7rc; TTt; 

Thus for Eq. [9] we have 

CUr 4 oc ADn 6 £- 4 , 

and for Eq. 0 we have 

1 f D 

CU£ -4 oc ADn 3 0 C 4 ■ 5a— / /dx. 

E Jo 

These two quartic divergent terms serve as the self-energy 
of the inter-media which is exactly of volume AD, and 
their independence of position a indicates they will not 
come into the Casimir force. Next, we also have 


which means CU£~ 3 = 0 for both Eqs. [T] and [9] This 
term is usually proportional to the surface area A and 
serves as the surface energy. The absence of surface diver¬ 
gence term is due to the cancellation of TE and TM con¬ 
tributions as also reported in Ref£. Next, we note, what¬ 
ever 1 1 £ 7 is, we have S(£) • lit; 1 = 0 and whatever Z 2 £ 2 is, 
we have c^S(£) ■ Z 2 £ 2 = 0. Thus C_ 2 £ -2 = C_i£ _1 = 0. 
Therefore, the cutoff-dependent terms in Eqs. [7] and [9] all 
vanish in the Casimir force. 

Now we turn to the contribution from Eq. [HI which has 
been expressed as Eq. [16] Since 9| has a good property— 
all polynomial terms under it vanish when £ —> 0, for 
positive plural s = 2?’, we only need to consider logarithm 
terms in Li s to check the cutoff property. We have 


»f(!) l0g « + S? (D^) l0g{ 



(23) 


After performing 9|, we see, term i = 1 contributes to 
C_ 2 £ -2 and term i = 2 contributes to Clo log£. All other 
logarithm terms under d | vanish when £ —> 0. We see 
C -2 oc jj, Clo oc -pj, and thus both of them are indepen¬ 
dent of position a. However, in half-space model there 
is no right cavity, those logarithm terms will appear to 
be a-dependent and must be removed manually. Unfor¬ 
tunately, this is an unresolved problem yet, as we kown. 

This completes our proof that for Taylor basis of any 
order d , the Casimir force for plates within an inhomoge¬ 
neous medium is cutoff-independent. Therefore, for any 
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smooth inhomogeneity, which is a superposition of Tay¬ 
lor bases, the Casimir force will have the inherited cutoff 
independence. 


V. APPLICATIONS 


The Taylor expansion of inhomogeneity profile, is use¬ 
ful for cutoff analysis, but will result in a series of con¬ 
stant terms Cq. This is inconvenient for calculation of 
the free Casimir energy and the Casimir force. Fortu¬ 
nately, for some cases, we are able to do the calculation 
without the Taylor expansion. One example is the profile 
given in Eq. [18] with 77 > 0, A = 0. 

The cutoff independence of Eh,5E\ can be analyzed 
exactly in the same way as above, while the 5E2 now is 
described by Eq. [20] Similarly, the logarithm terms 


(1 loglog? 

oc (1 — e~ vD )£ k log£, 


is independent of position a. Therefore, we have again 
the cutoff independence of Casimir force, as expected. 

To obtain the explicit expression of the Casimir force, 
we only need to calculate the Co term, which comes from 
£ 3 in Li or </>. Now we have 


1 


F 0 C = — (1-e-n, iff = m , 
rja rj(D — a) 


{e~ va - e~ vD ). 


Cq contribution from Eh is 


AlS(-2) 

47T 

AHctt 2 


2C(—3) 


2C( 3) . 7T 3 

3! ano ' 3! v (D — a)no 


(-- 


(24) 

) 3 


1 


1 


720n o La 3 (D - a ) 3 ' 


(25) 


Here we have recovered He in the end. Co contribution 
from 6 Ei is 


5a 


AHctt 2 

720no 


F 0 C (-) 3 + 

a D — a 


(26) 


During the calculation of this term, we have seen the 
operator dc ■ £• introduced by including the exponent in 
the regularization in Eq. [3] does not change the observ¬ 
able value (compared with the one without including the 
exponent in the regularization), as expected. Cq contri¬ 
bution from 5E2 is 


1 180 

72(Jno a 7 r 


/3 3 I ^( 1 ) - log + t 1 ) + /3 3 (log i/3 - ijj{/3)) - ,9/12 


11 


2 n n 


+ £ -A n, 


where /? = — ^ and C — > 1Z represents the right counter¬ 
part where we should make a replacement (1 — e ~ r]a ) —> 
[e~ va — e~ vD ) and replacement a —> {D — a) for other 
a’s. S is the symbol for imaginary part of an expression. 

If we make a transform rj —>■ —*5-Jy, A —> iA and make 
use of the real part of Eas. [Ihlandl^Dl we can have some 
insights for profiles f(x) = cos (^x + A), where 0 < 
A < 27 t and b > 0. We have 

F ™,b = Nn A - sin (2 tt|^ + A)] 

1 1 
m — ba/2 tt m + ba/2n 


and 


r ~ 5aAn?, 1 r . A . . ba 

SE2 = 27t 4tt Sln A ~ sm(27r 2 D + A)] 




0 (e“ 




! ba 
n ° , 1 . — 

’ ’ 2tt‘ 


Here we choose 5 = 2 and A = 0. For 0 < a < D we 
have 0 < ba/2D < 1 so the Lerch function is well defined. 
Together with another profile /( x) = 1 and some numeric 
factors, we can recover the result of Ref. 12 - following the 
procedures above. 

Back to the exponentially decaying profile, we evaluate 
the influence of inhomogeneity on the total Casimir force. 
We focus on the left cavity part. When a —> oo, 5ft[log i/3— 
Tp{f3)\ vanishes, but ^5(1) — log ^ + A j i s nonzero and 

depends on materials properties no, rj. This means Ea.l27l 
contains a part of energy that is not free and thus has no 
influence to Casimir force. We let D —> oo to remove the 
right cavity so that we can focus only on the left cavity 
(two-plate interaction). The force contributions are 


F h = - 


5Fi = 5a 


AHctt 2 1 
240?r o a 1 ’ 

F 0 C + l(F 0 C - /(a)) 


AHctt 2 1 

240n o a 1 3^’° 


(28) 

(29) 


6F 2 = 5a 

60 

TT 


AHctt 2 1 f 60 


240?ro a 4 ] 7 r 


(e _r?a - 1 )S 


/ 3 3 (1 — P4>'(/3)) + 


Vae-^Z /? 3 (<M 1 ) - log ^ 


P 


+ W (log i/3 — ip(P)) ~ Y 2 


(30) 


Eq. [28]is the well-known Casimir force between two plates 
within homogeneous media. If we treat the medium be¬ 
tween two plates as homogeneous and use the average 
refraction index, we can get an approximation of the 
Casimir force 


F h = 


AHctt 2 


1 


[1 -SaFf]. 


(27) 


240?ro a 4 


(31) 
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FIG. 2. The homogeneous approximation and the exact 
Casimir pressure including the first-order correction between 
plates within an inhomogeneous medium. Solid curves are ap¬ 
proximations and dashed curves are exact pressures in unit of 
Pascal (minus means attractive). All x axis is the position a. 
The parameters are rj = 10 7 /m,no = 1.5 for all, while 5a = 
0.3, 0.1, 0.01 for red (long-dashed), blue (medium-dashed) and 
cyan (short-dashed) curves respectively. Red solid, blue solid, 
cyan solid and cyan short-dashed curves appear to coincide. 
Insets are ratios of 5F/Fh . Top inset: 5a = 0.1, no = 1.5 and 
r/ = 2 x 10 7 ,5 x 10 6 ,2 x 10 6 ,5 x 10 5 /m for red solid, green 
short-dashed, blue medium-dashed, cyan long-dashed curves, 
respectively. Bottom inset: rj = 2 x 10 7 /m, no = 1.5 and 
5a = 0.3,0.2,0.1,0.01 for red solid, green short-dashed, blue 
medium-dashed, cyan long-dashed curves, respectively. 


This is exactly the combination of Eq. [28] and the first 
part ofEq.ED The rest (second part) of Eg. IHTTI together 
with Eq. [3D1 is written as 5F. This term is easy to un¬ 
derstand. It reflects the change of the average refraction 
index when the plate is shifted. 

The relation between Fh and 5F is given in Fig. [2] 
We should emphasize that, according to Eqs. [5] and HOI 
we have 5a^2 x Vj < 2 da. All 5a < 1 is permitted 
within perturbation theory (though the first-order cor¬ 
rection might not be enough). Our simulation results 
clearly show, when ?? > 10 , 5a > 0.2, the correction even 
dominates over the homogeneous approximation and flips 
the sign in some range. When 5a = 0.1, the first- 
order corrections also have a relative magnitude of peak 
50%, thus can not be ignored. We should take seriously 
this inhomogeneity-induced repulsion, since this might 
indicate alternatively the first-order perturbation is too 
rough for such an intense inhomogeneity. To investi¬ 
gate whether inhomogeneity can induce repulsive Casimir 


force and be used to control the Casimir force, we need 
further studies. On the other hand, when 5a -C 0.01, 
or rj -C 5 x 10 5 /m, in the range of 0~1 um where the 
Casimir force is measurable, the correction is well below 
1% and can be omitted. Actually, in our minds, most 
experiments in natural condition are belong to this case. 

VI. DISCUSSIONS 

First, we note the linearly inhomogeneous case f{x) = 
x/D. From Eq. [13] we know 5E2 = 0, thus there is no 
logarithm term. This seems to tell us that, removing the 
“bulk contribution” is enough in this case to retrieve fi¬ 
nite Casimir force in the half-space model. This actually 
is not true. A second-order correction ( Xj\( x /D) 2 \xj) 
would immediately introduce in logarithm terms, let 
alone higher order corrections. 

We thus stress that our result is valid within the 
first-order perturbation theory, because, without a te¬ 
dious calculation of higher order perturbation or any 
other rigorous demonstration, one can never claim there 
wouldn’t be any omitted term that is actually both in¬ 
dependent and cutoff-dependent. One example is, 5E ~ 
£ -4 + alog£ + /(a). When £ —> 0 + , to calculate per¬ 
turbative energy for the first-order perturbation, alog£ 
may be possibly ignored as a trivial term. Since £ -4 is 
independent of a, it will disappear in the expression of 
Casimir force. Then everything looks fine and no cut¬ 
off dependence is found in such an approach. However, 
what matters the most to the Casimir force is the deriva¬ 
tive of energy with respect to a, and if we know the exact 
perturbative energy somehow, then we see the previously 
trivial term alog£ now arises to be cutoff-dependent. 

At last, if somehow we know higher order corrections 
would not produce both a and cutoff dependent terms, 
which is a reasonable expect, the magnitudes of first- 
order corrections would be quite reliable then. 

VII. CONCLUSIONS 

With mode-summation technique and first-order per¬ 
turbation theory, we have expressed the regularized 
Casimir energy for inhomogeneous case with cylinder ker¬ 
nel coefficients, as Eq. [23] Like other unresolved Casimir 
divergences, we found the presence of the logarithmi¬ 
cally cutoff-dependent term (see Eq. [33] and the subse¬ 
quent analysis). Our results also show there is a term 
of quadratic cutoff dependence in the Casimir energy. In 
the piston model such terms are independent of position 
a and thus vanish in the force, while in the half-space 
model such terms are dependent of a and thus remain 
in the force. Consequently, we must introduce additional 
regularizations to remove them in the half-space model, 
though, it is not clear how to do it. 

Based on the piston model, our results have shown, for 
any smoothly inhomogeneous profile, the Casimir force 
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is always cutoff-independent in the first-order perturba¬ 
tion. For some other profiles that are not smooth, it 
seems one can still get cutoff-independent result, though 
we can not give a rigorous proof that it’s always the case 
yet. Our result supports the method in RefJ^ to omit di¬ 
verging terms when simulating the Casimir force within 
inhomogeneous media numerically. 

We have also calculated the first-order corrections to 
both free Casimir energy and Casimir force for exponen¬ 
tially decaying profile. Surprisingly, comparing with the 
homogeneous analogue where the average refraction in¬ 
dex between two plates is used, we found the correction 
to Casimir force can be even larger than the predicted 


value of the homogeneous analogue, and flips the sign of 
the force, though we note the first-order correction might 
not be accurate enough. All of these results may be useful 
as a reference for future relevant theoretical calculations 
and experimental studies. 
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